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Abstract. - A numerically exact dynamical quantum Monte Carlo approach has been de- 
veloped and applied to transport through a double barrier in a Luttinger liquid with arbitrary 
transmission. For strong transmission, we find broad Fabry-Perot Coulomb blockade peaks, 
with a lineshape parametrized by a single parameter, but at sufficiently low temperatures, non- 
Lorentzian universal lineshapes characteristic of coherent resonant tunneling emerge, even for 
strong interactions. For weak transmission, our data supports the recently proposed correlated 
sequential tunneling picture and is consistent with experimental results on intrinsic nanotube 
dots. 



Resonant tunneling in a (non-chiral) Luttinger liquid (LL) was studied more than a decade 
ago [1-3], but has recently attracted widespread attention by theorists again [4-10]. This 
interest is primarily caused by new experimental realizations of double-barrier setups in inter- 
acting ID quantum wires presumably described by LL theory, using semiconductor quantum 
wires [11] or nanotubes [12]. The experiments of Ref. [12] have been interpreted in terms of 
a "correlated sequential tunneling" (CST) mechanism [7,12], since the standard picture of 
uncorrelated sequential tunneling (UST) in a LL [4] is inconsistent with the observed temper- 
ature dependence of the conductance peak height. However, as CST theory essentially relies 
on a master equation [13], it is clearly of interest to check it against exact results. In this 
context also other nanotube experiments are of interest, where for nearly transparent double- 
barrier, Fabry-Perot oscillations in the gate- voltage dependence of the conductance have been 
reported [14-17]. 

As the double-impurity problem in a LL is not integrable, exact solutions covering a wide 
parameter range are out of reach, and analytical progress has to rely on approximations. 
One line of reasoning considered very weak Coulomb interactions, i.e. Luttinger parameter g 
very close to one [8,9], where no CST processes were found. Furusaki has studied the UST 
regime, where the linewidth of the resonance peak has a linear temperature dependence, and 
the peak conductance Gp cx T~'^~^^^^ [4]. The CST mechanism could provide an additional 
transport channel that may dominate on resonance for strong interactions, and gives instead 
the observed [12] behavior, Gp oc T-3+2/9 [7]. The CST mechanis m is similar but different 
from conventional cotunneling, which can be the dominant transport channel away from the 
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resonance. We mention in passing that the exact 5=1/2 solution of a modified model in Ref. 
[10] neither contains the UST nor the CST regime and therefore cannot resolve this debate. 
Finally, at low temperatures, instead of sequential tunneling, coherent resonant tunneling is 
possible, characterized by non-Lorentzian universal lineshapes [1], 

G{N^,T)/G^^ fg{X), X = c|iVo-l/2|/ri-^ (1) 

with dimensionful constant c and Go = /h. Here A^o is the dot's average charge induced by 
external gates, with peak center, say, at iVo = 1/2. Obviously, the linewidth then scales as 
T^^3. The universal scaling function has the limiting behavior [1] 

1-/<,(X«1)«X2, /<,(X » 1) « (2) 

While for strong barriers, coherent resonant tunneling is only expected for g > 1/2 [4], for 
weak barriers (strong transmission), it is largely unclear to what extent this concept applies. 

Here we present numerically exact results obtained from a dynamical quantum Monte Carlo 
approach that provide detailed insight into the strong transmission regime and also sheds light 
on the controversy about mechanisms of resonant tunneling. This method was successfully 
used for the corresponding single-barrier case [18,19], and is generalized here to resonant 
tunneling. We focus on the linear conductance for spinless electrons and symmetric barriers, 
since additional simulations show that neither spin/flavor degeneracy nor weak asymmetry 
have a qualitative effect on our results. 

We consider a LL containing two impurities of strength Vb at x = -iid/2, thereby forming 
a quantum dot with LL leads. The single-particle level spacing is Es — nvpjd, the charging 
energy is £^c = Es/ ■ In terms of the standard boson field (j){x) and its conjugate momentum 
H(x), the Hamiltonian is [1] 



I dx I h2 + ^ (9^0)2 1 + Vb ^ cos[V4^,/.(pd/2, t) + eVt + pnNo] , (3) 

where vp is the Fermi velocity, V the applied bias voltage, and the current through the dot is 

I^GoV+^{dtHx,t)), (4) 

where x is arbitrary and t oo. For g = 1, refermionization yields the exact conductance 
G = dl/dV. With bandwidth D, the dimensionless linewidth 

^^8A(I+i)' ^ = -VA (5) 
and the derivative of the Fermi function, -df/dE = l/lAT cosh^{E/2T)], we obtain 



Go 7-00 dE cos2(7r[7Vo + E/Es]) + ' 

where h = ks = 1- For strong barriers, this leads to the usual Breit-Wigner lineshape with 
linewidth wEs/n. Note that the infinite-barrier limit is reached already for A = 2, where the 
associated phase shift is in the unitary limit. Equation J^Jl holds for arbitrary barrier height 
Vo, including strong transmission (Vb — > 0), and allows to firmly establish the validity of our 
numerical scheme. 
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Next we outline our path-integral Monte Carlo (PIMC) approach to the linear conductance 
for arbitrary g. While PIMC is conventionally used to evaluate imaginary-time path integrals, 
conductance calculations need dynamical information. We first tried to use various schemes 
to analytically continue imaginary-time PIMC data, but the results were not reliable. This 
reflects a well-known difficulty related to the numerically ill-posed nature of the analytic 
continuation [19]. We therefore proceed directly within a (Keldysh) real-time formalism. 
Although real-time PIMC has to deal with the sign problem, our formulation avoids much of it 
by mapping the problem to an equivalent Coulomb gas description. In this representation, the 
sign problem is rather weak and permits numerically exact simulations for the full parameter 
regime of interest. 

Consider the discretized Keldysh contour running from time t ~ Q to t,nax and back 
to zero. We keep t,nax finite and define a discrete time spacing A = imax/-P with Trotter 
number P. At time tj = [j — 1)A, fields (t)j{x) (0^(x)) live on the forward (backward) 
branch. Next we switch to a Coulomb gas picture by expanding the impurity propagator 
for sufficiently small AVb [19]. Following Ref. [18], we use a Coulomb gas expansion valid 
up to order (AVb)^. Introducing "quantum" charges ^jp = 0, ±1,±2 and "quasiclassical" 
charges rjjp ~ 0,±1/2,±1 for each time (j — 1, . . . , P) and impurity index (j) = ±), where 
77 ± ^/2 must be integer, it is sufficient to keep |77 + ^/2| + [77 — ^/21 < 2 within this order 
of accuracy. Only configurations subject to electroneutrality, X^jp ^jp ~ contribute to the 
partition function. Moreover, it turns out that the quasiclassical 77 charges can be summed over 
analytically [18,19]. With z defined in Eq. Q below, this leads to effective Greens functions 
K{^,z), with the entries i^(0, z) = 1 - 2(At/o/2)2(l - cos 2), K{±l,z) = ±AVb sin(2/2), and 
K{±2,z) — (AVb/2)^(l — cosz). Under the Coulomb gas expansion, we can now integrate out 
all boson fields 4>i{x) and 4''i{x), since they appear only quadratically in the action. Thereby 
we arrive at an effective action governing the dynamics of the Coulomb gas charges {^}. The 
result can be put into the language of dissipative quantum mechanics [20] by defining spectral 
densities J±{ijj) — 7r(7ti;[l ± cos{t: lo / Ec)]e~'^ ^ , with associated correlation function 



This implies (i) the action contribution <i>'[{C}] = Spp' I]j>fe Cjp['5'+ jfc + pp'S*- jfc]Cfep'; and 
(ii) the z's entering if (^, z) are given by 




00 



duj J±(w) cosh[w(l/2r- ii)] - cosh[cj/2r] 
TT sinh[ti;/2T] 




(7) 



j>k p' 



with S±^jk + iR±jk — L±{[j — k]A); for the diagonal elements, see Ref. [21]. 
Collecting results, the conductance is obtained as 

G(iVo,T)/Go-l - lim dtlsit), 



(8) 



where the function Ib (t) can be computed from 



lB{tk) = Z-i^/fe[U}]cos TTiVo^pOp exp(-$'[{e}])n^(Op.^.p)- 




(9) 



The normalization Z is the path sum for Ik 



1, and 



j'p' ip j<k,p 



■+.kj 



9zK{£_jp, Zjp) 
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Fig. 1 - Linear conductance versus iVo for g = l,Es/D = ■k/2,T/D = 0.025, with Vo/D = 0.05 
(PIMC: circles; Eq. ||SJ: solid curve), Vo/D = 0.1 (squares; dashed curve), and Vo/D = 0.2 (diamonds; 
dotted curve). 



This is a real- valued quantity, as are all the other quantities appearing in Eq. Remarkably, 
although we are dealing with a real-time sign problem, it effectively looks just like a fermion 
one. The sign problem arises because the integrand in Eq. ® can be negative, leading to 
interference effects and small signal-to-noise ratio at sufficiently long real times imax- 

Monte Carlo trajectories for the quantum fluctuations {^} are then generated by taking the 
absolute value of the integrand in Eq. @ as weight. The average sign in the data reported here 
was always larger than « 0.001, which still allows to run stable simulations, albeit sometimes 
at the expense of very long CPU times. One can obtain the whole function I sit) in one MC 
run. For sufficiently long times, this function has a well-defined linear slope which determines 
the conductance via Eq. ©. Typically, Trotter convergence was reached for discretizations 
AVq < 0.1. On a 2 GHz Xeon processor our code performs at an average speed of about 10^ 
samples per hour (for P = 40). Several 10^ MC samples have to be accumulated to obtain 
Isit) for a given parameter set in order to ensure good statistics. Error bars then refer to 
both standard stochastic MC errors and to uncertainties from fitting the long-time behavior 
by a linear slope. The validity and accuracy of our scheme has been established by checking 
numerical data against the exact 5 = 1 solution, Eq. see Fig.^ The comparison highlights 
the power of our approach. We then move on to interacting electrons, focusing on g = 0.3 
and g = 0.6. 

Let us start with the case of strong transmission, taking V^/D = 0.05. Remarkably, for 
T/D > 0.01, PIMC data are quantitatively described by the g — I lineshape © provided 
w — Wg{T) is treated as a fit parameter. The corresponding values of w are shown in Fig. |21 
and reveal power-law behavior, Wg{T) oc T^f, with g-dependent exponent ag. For very 
weak interactions, a related behavior has been discussed before [9]. Here we find numerical 
evidence for a power-law temperature dependence of the linewidth for strong interactions. For 
the two values of g studied, we obtain ao.3 ~ 0-84 and ao.e ~ 0.72. Therefore the strong- 
transmission peaks become narrower and narrower as T is lowered. Note that each data 
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Fig. 2 - Fit parameter 'Wg{T) for conductance lineshape @ as obtained from QMC data a.t g = 0.3 
(diamonds) and g — 0.6 (squares), for Es/D — 7r/2 and Vo/D — 0.05. Dotted and dashed lines are 
guides to the eye only. 

point shown in Fig. |21 actually has been obtained from QMC data for the full conductance 
lineshape. The lineshape © closely resembles experimental results for strong-transmission 
(Fabry-Perot) oscillations in nanotubes [14-17]. We therefore identify this region of weak 
barriers and not too low temperature as the Fabry-Perot regime. For g < 1, such Fabry- 
Perot oscillations include remnants of Coulomb blockade effects, which are partly washed out 
due to pronounced quantum fluctuations present at strong transmission. Nevertheless, these 
effects are responsible for the narrowing of the resonance peak as temperature is lowered. 
We mention in passing that no narrow conductance dips were observed such as the ones seen 
experimentally in Ref. [14]. Such dips are probably related to special impurity scattering 
processes not contained in the model iQ. 

At lower temperatures, deviations from the Fabry-Perot lineshape (0 can be seen. How- 
ever, then our data can be collapsed onto the universal scaling curve Q), see Fig.^lfor g = 0.3. 
Very similar results were found for g — 0.6 as well. For small and large X — c|iVo — 1/2|/T^~^, 
the results in Fig. Oare consistent with the respective analytical prediction Therefore 
these universal lineshapes can be identified as coherent resonant tunneling peaks. Although 
for strong barriers, resonant tunneling exists only for g > 1/2 [1], we observe a perfect reso- 
nance peak a,t g = 0.3. This is in accordance with renormalization group arguments for weak 
impurities for g > 1/4 [1], and shows that the picture of coherent resonant tunneling in a 
Luttinger liquid is actually very robust. Only for very weak barriers and sufficiently high T, 
the Fabry-Perot regime replaces the universal region of resonant tunneling. 

Next we discuss stronger barriers, Vo/D = 0.2, where one expects sequential tunneling at 
not too low temperatures. In Fig. 01 data for the temperature dependence of the conductance 
peak height, Gp{T), is shown for g = 0.6. At low temperatures, Gp approaches the perfect 
quantum conductance Gq expected in the resonant tunneling regime. With increasing tem- 
perature, the conductance goes through a minimum. For intermediate temperatures, but still 
well below Es, sequential tunneling then starts to dominate, and Gp(T) increases in a power- 
law fashion, Gp{T) oc T"^ . Remarkably, our data are consistent with the CST scaling [7], 
77 = 2/(7 — 3 = 1/3, but not with the corresponding UST law [4],77 = l/5 — 2 = — 1/3, which 
would even have a different sign. We stress that other exponents, e.g. 77 = 2/(7 — 2 = 4/3 as 
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Fig. 3 - Low-temperature QMC data for various No,T at Es/D = tt/2,Vo/D = 0.05, plotted 
according to Eq. Q for g = 0.3. Different symbols refer to different temperatures. 



expected for the effective strong single-impurity problem arising at T 3> i?s are incompati- 
ble with this numerically exact result. Therefore our simulation data suggests that additional 
mechanisms beyond conventional sequential tunneling can indeed be crucial. As a side remark, 
we mention that Fig.^shows similar high-temperature (T « Es) features as the experimental 
data of Ref. [12]. The power-law upturn for Gp{T) towards the resonant tunneling limit, see 
Fig. 0] at T/Eg < 0.1, has not yet been observed experimentally [12], probably due to a larger 




Fig. 4 - Temperature dependence of the peak conductance Gp for g — 0.6, Vo/D = 0.2, and 
Es/D = n/20. The dashed line reflects a T"" * power law, the solid line a T^^^ law. 



S. HUGLE and R. EGGER: RESONANT TUNNELING IN A LUTTINGER LIQUID 



7 



barrier height Vq in the samples. 

To conclude, we have developed and applied a numerically exact and well-controlled real- 
time Monte Carlo approach to the computation of the resonant tunneling conductance in a 
Luttinger liquid. For weak barriers and not too low temperature, we identify a Fabry-Perot 
regime, where the conductance peak has a lineshape given by Eq. JHl with a temperature- 
dependent linewidth parameter Wg{T). Within the range of applicability, Wg(T) exhibits 
g-dependent power-law scaling. At sufficiently low temperatures, for g > 1/4, we then find 
a crossover into the universal coherent resonant tunneling regime. At higher temperatures, 
we find numerical evidence in support of the recently suggested correlated sequential tunnel- 
ing picture [7]. Finally, we note that this algorithm also allows to study out-of-equilibrium 
transport, quantum shot noise, and transport in the presence of a Kondo-type dot. 

* * * 
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